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Abstract 

The coincidence detection efficiency of the TUNL-ITEP apparatus designed for 
measuring half-hfe times of two-neutrino double-beta {2u(5(5) decay transitions to 
excited final states in daughter nuclei has been measured with a factor of 2.4 im- 
proved accuracy. In addition, the previous measuring time of 455 days for the study 
of the ^'^'^Mo 2v(3l3 decay to the first excited 0]'' state in ^'^'^Ru has been increased 
by 450 days, and a new result (combined with the previous measurement obtained 
with the same apparatus) for this transition is presented: Tx/2= [5.5lQ'g (stat) ± 
0.3 (syst)] X lO^'' y. Measured 2v(3() decay half-life times to excited states can be 
used to test the rehabihty of nuclear matrix element calculations needed for deter- 
mining the effective neutrino mass from zero-neutrino double-beta decay data. We 
also present new limits for transitions to higher excited states in ^'^'^Ru which, if 
improved, may be of interest for more exotic conjectures, like a bosonic component 
to neutrino statistics. 
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1 INTRODUCTION 



The half-life time of the zero-neutrino double-beta {0i>(3j3) decay is related to 
the neutrino mass by the equation 

[T?/2r' = GoAM'/ - igA/gvfM'Jry <m,>' • (1) 

Here, Goiy is the two-body phase-space factor which includes the coupling con- 
stant, <mi, > is the effective Majorana electron neutrino mass, and M^^ and 
M.o'rp are the Fermi and Gamow- Teller nuclear matrix elements, respectively. 
Therefore, extracting the neutrino mass from the measured half-life time re- 
quires the knowledge of the Oz//3/3 nuclear matrix elements. Information on 
the nuclear matrix elements is currently obtained from quasi-particle random 
phase approximation (QRPA) and nuclear shell-model calculations pLj. The 
parameters in these calculations, in particular the strength of the particle- 
particle interaction gpp, in turn are adjusted using measured half-life times of 
single-beta, and more recently, two-neutrino double-beta (2z//3/5) decay data. 
Because the behavior of gpp is different for the transition to the ground state 
compared to transitions to excited states in the daughter nucleus, measuring 
2z//3/5 decay to excited final states provides an excellent test of the theoret- 
ical calculations [2]. Furthermore, a very interesting possibility is the search 
for the OuPP transition to the 0^^ excited state of the daughter nucleus. This 
transition provides a clear-cut signature: in addition to the two electrons with 
a fixed total energy, there are two gamma rays whose energies are strictly fixed 
as well. In a hypothetical experiment detecting all decay products with high 
efficiency and sufficient energy and spatial resolution, the background can be 
reduced to nearly zero. It is quite possible this idea will be used in future large 
scale 0up(3 searches such as MAJORANA [3], [1], GERDA 0, [S], CUORE 
[7], [8], or SuperNEMO [9],[l0]. In addition, in Ref. [H] it was stated that 
the detection of the excited state transition will provide the opportunity of 
distinguishing between different OuPP mechanisms, such as light and heavy 
Majorana neutrino exchange or trilinear R-parity breaking mechanisms. 

Recently it was speculated [12] that neutrinos may violate the Pauli exclusion 
principle (PEP) and therefore, at least partly obey Bose-Einstein statistics 
(see also [I3]). As a consequence, neutrinos could form a Bose condensate 
which may account for parts or even all of the dark matter in the universe. As 
discussed in [14j the possible violation of the PEP has interesting consequences 
for 2i'PP decay. It not only modifies the energy and angular distributions of the 
emitted electrons, but it also strongly affects the 2i//?/3 decay rates to ground 
and excited states in daughter nuclei. Following [T3], the half-life time ratios 
for transitions to excited 2+ states and the 0^ ground state are by far the most 
sensitive way of obtaining bounds on the conjecture of a substantial bosonic 
component to neutrino statistics [12]. As a result, information on the decay 
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rates to excited states is needed to test this new and potentially far-reaching 
hypothesis. 

Presently, 2i//5/5 decay has been observed for ten nuclei [15], but data for the 
2upp decay to excited final states exist only for i°°Mo [l6]-[20] and ^^°Nd [21] ■ 

Recently, a new determination of the half-life time of ^°'^Mo for the transition 
to the first excited 0^^ state in ^°°Ru was reported by Hornish et al. [19]. In 
addition, limits to higher excited states were also given in [19j. These mea- 
surements were performed on ground level at TUNL using the coincidence 
technique. In the case of ^°°Mo, the excited state of the daughter nucleus 
^°°Ru de-excites via the decay sequence Qi^2i^ 0^^, emitting two 7 rays 
in coincidence with energies of E^i=590.8 keV and E^2=539.5 keV. 

The TUNL-ITEP (3(5 decay setup consists of two cylindrical high-purity ger- 
manium (HPGe) detectors with the (3(3 decay source sandwiched between 
them. The HPGe detectors are surrounded by a Nal(Tl) annulus which in 
turn is placed inside of a lead-brick enclosure. The combination of active and 
passive shielding suppresses the coincidence background to such an extent 
that successful measurements on ground level are feasible, provided (3(3 decay 
source masses of about 1 kg are available and half-life times do not exceed 
10^^ years. The (3(3 decay source used at TUNL was a metallic disk with mass 
of 1.05 kg, enriched to 98.4% in ^°°Mo. After counting for 455 days, 22 coinci- 
dence events were detected with a continuous background estimated to be 2.5 
events. Using the measured efficiency of the detector system, a half-life time 
of Ti/2=[6.0il:?(stat)±0.6(syst)]xl020 y was reported in [19]. 

Here we report an improved measurement of the coincidence efficiency of the 
TUNL-ITEP (3(3 decay apparatus, and new data taken for ^''^Mo during 450 
days of counting. 



2 COINCIDENCE-EFFICIENCY MEASUREMENT 



We closely followed the procedure described in detail in [TH]. The radioactive 
source used to measure the efficiency was ^°^Rh, which was produced at TUNL 
via the reaction ^°^Ru(p,n)^°^Rh initiated with 5 MeV incident protons on a 
natural ruthenium target (31.6% ^''^Ru). This particular source was chosen 
for several reasons. First, the O^*" ^ 2^^ ^ 0^^ decay sequence not only results 
in 7 rays in coincidence, but also mimics the decay scheme of the excited 
final 05^ state of ^°°Mo. Therefore, corrections for the angular distribution are 
unnecessary. The associated 7 ray-energies, 468.6 keV and 475.1 keV, are also 
close to those of ^°°Ru (539.5 keV and 590.8 keV), the daughter nucleus of 
^°°Mo. Finally, the 205.1 ± 1.9 day half-life of ^"^Rh is long enough to allow 
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shorter lived contaminants to decay away prior to measurement. The activity 
of the source was determined by comparing the intensity of the 475.1 keV 7- 
ray peak to that of the 661.6 keV 7-ray peak from a cahbrated ^^''Cs source. 

The HPGe detectors of [19j were 8.8 cm in diameter and 5.0 cm in thickness. 
Because the ^°°Mo source used for the measurement of the half-hfe time in [19] 
and in our work described in Sec. 3 was in form of a 10.6 cm diameter and 1.1 
cm thick disk, the radial dependence of the 468.6 keV - 475.1 keV coincidence 
efficiency across the detector faces was crucial. We sandwiched our 1 x 1 x 0.5 
mm^ ruthenium/rhodium source (in the following referred to as ^°^Rh source) 
between ten disks (five on each side) of 10 cm diameter and 0.1 cm thickness 
of natural molybdenum to ensure that the overall thickness was the same as 
in the original 2i>(3f3 decay experiment on ^°'^Mo. This procedure yields the 
same average 7-ray attenuation provided the energy difference between the 
7-ray pairs (468.6 keV - 475.1 keV versus 539.5 keV - 590.8 keV) is properly 
taken into account. Furthermore, it enables the ^°^Rh source to be positioned 
not only at different radii, but also at different locations along the z axis, i.e., 
along the common axis of the two HPGe detectors. 

The previous coincidence efficiency measurement [12] was made along the 
radius r on the front faces of the HPGe detectors at the four locations r=0, 
2, 4 and 5 cm, i.e., in only one direction, assuming cylindrical symmetry of 
the detectors. The current measurements were performed across the entire 
diameter and in two perpendicular directions, i.e., horizontally and vertically. 
Near the center we used 0.5 cm increments to better determine the efficiency in 
that region. A total of 25 data points were obtained in this manner between 
-4.5 cm < r < 4.5 cm with the ^°^Rh source sandwiched between two 0.5 
cm-thick molybdenum disks. The measuring time for each data point varied 
between 12 and 48 hours, always obtaining at least 8000 coincidence events. 
At each location of the ^°^Rh source the number of 468.6 keV - 475.1 keV 
coincidences was determined from the two-dimensional data area of pulse 
height in HPGe detector 1 versus pulse height in HPGe detector 2 (see Fig. 

The data were analyzed using the ROOT C/C++ Interpreter. The program 
used gates and projections (see Fig. [2]) along the axes to set close gates around 
the coincidence peaks and, after background subtraction, obtained the yields 
in each peak. The time normalized yields were then corrected for dead time 
(<3%) and normalized to the activity of the ^°^Rh source, taking into account 
its half-life. We measured this half-life independently and determined it to be 
203.3 ± 2.8 days. Combined with an existing world average of 206.7 ± 2.6 
days from [22]-[23], we obtained a half life of 205.1 ± 1.9 days. In order to 
obtain the efficiency for the ^°°Mo 2z//3/3 decay to the 0^ state of ^°°Ru the 
yields had to be corrected not only for the difference in attenuation between 
the ^°^Rh and ^°°Ru 7-ray lines of interest, but also for the energy dependent 
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detection efficiency of our two identical HPGe detectors. For that purpose, the 
relative detection efficiency was measured with a ^^^Eu source. 

As has been expected, the radial dependence of the coincidence efficiency was 
almost identical for the horizontal and vertical scans, within error associated 
with the exact location of the ^°^Rh source. Figure [3] shows the average of the 
horizontal and vertical efficiencies as a function of r. The error bars represent 
the statistical uncertainty of the data. Here, the effect of the position uncer- 
tainty of the ^°^Rh source is added in quadrature to the statistical uncertainty. 
The coincidence efficiency is rather small, close to 0.4% in the -1 cm < r < 1 
cm range and then dropping smoothly to about 0.1% at r=4: cm. As a result, 
only about 3% of our ^°°Mo disk provides 2v(3(3 decay related 7 rays with the 
maximum detection efficiency. However, about 30% of the disk provides 7 rays 
with an efficiency of 50% of the maximum efficiency obtained at r=0 cm. 

The curve through the data points in Fig. [3] is a least-square fit using the 
functional form 



where a, 6, and c are free parameters and r^O refers to the center on the front 
face of the HPGe detectors. A careful inspection reveals a slight asymmetry 
in the coincidence efficiency, providing slightly larger values for r< -2 cm. 
The lower (dashed) and upper (dashed-dotted) curves shown in Fig. [3] repre- 
sent our systematic uncertainty of 5.1%, including the ±3% scale uncertainty 
associated with the "^'^'^Cs 7-ray source. 

The radial contribution to the coincidence efficiency was calculated using 
/TT J rar 



where e^-y(r) is the best fit obtained from an asymmetric fit to the data. This 
value is then corrected by the ^-dependence of the coincidence efficiency mea- 
sured in the present work and confirmed by Monte-Carlo simulation, where 
z is the distance from the face of the detector. There is a 10% decrease in 
efficiency at the center of the 1 cm thick molybdenum disk compared to the 
front and back faces. The resulting efficiency value is etot=(0.182±0.009)%, 
where the total uncertainty of 5.1% is due to the contributions listed in Ta- 
ble [H They include an estimated 3% uncertainty due to the slightly irregular 
shape of our ^°°Mo disk. The previous measurements described in [T9] provided 
etof=(0.219±0.022)%. The associated data are shown in Fig. H] in comparison 
with a fit to the present data (solid curve). The main reason for the difference 
between the previous and the present value is due to an improved determina- 
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tion of the dependence of the efficiency. 



3 NEW DATA FOR DOUBLE-BETA DECAY OF ^o^Mo TO THE 
FIRST EXCITED 0+ STATE OF io°Ru 



The improvement by a factor of 2.4 in the determination of the coincidence 
efficiency for the 2z//3/3 decay apparatus of [19] justifies reducing the statistical 
uncertainty of Ti/2=[6.0ti;?(stat) ±0.6(syst)] x lO^o y obtained in [L9\ for the 
half-hfe time of ^°°Mo for the transition to the first excited Of state in ^°°Ru. 
Only in the very recent result of the NEMO group of Ti/2=[5.7tJ;9 (stat) ± 
0.8 (syst)]xlO^° y [20] is the statistical uncertainty appreciably smaller and 
approaches the systematic uncertainty of the experimental data. However, in 
order for 21/(3(3 decay half-life time measurements to the 0^ state to play a 
decisive role in testing nuclear matrix element calculations, statistical and 
systematic accuracies each approaching the ~5% level are needed for ^°°Mo. 

In an attempt to approach this level of accuracy also for the statistical uncer- 
tainty, we report in this section on a new result of our ongoing measurement of 
Ti/2 for ^'^°Mo to the Of state in ^°°Ru. The part of the level scheme of ^°°Ru 
which is of interest to the present work is given in Fig. [51 Because the approach 
and analysis procedures are the same as those described in [19] , only the essen- 
tial facts are given here. After 450 days of data taking we recorded 16.0 ± 4.4 
net coincidence events after subtraction of 0.45 ± 0.16 background events per 
keV. This result can be compared to the 19.5 ± 4.7 net events reported in [19] 
for 455 days of data taking. Figure [6] shows the energy spectra for the 539.5 
keV and 590.8 keV transitions obtained in the present work, while Fig.[7]shows 
our spectra combined with those of [19], resulting in a total of 35.5 ± 6.4 coin- 
cidence events. Using our newly determined coincidence efficiency discussed in 
Sec. 2 we obtain Ti/2=[6.0fli (stat) ± 0.3 (syst)]xl02° y for the present data 
set. We have no indication that the efficiency of the TUNL 21/(3(3 decay appa- 
ratus has changed between the data of [19] and the present data were accumu- 
lated. Therefore, we used our newly determined coincidence efficiency for the 
data of [19] and obtained Ti/2=[5.0;o:9 (stat) ± 0.3 (syst)]xl020 y. Combining 
this result with our present result gives Ti/2=[5.5lJ;8 (stat) ± 0.3 (syst)]xlO^° 
y. This value is in agreement with the recommended value of Barabash [25] 
for the decay of i°°Mo to the 0^ state in i°°Ru of Ti/2=6.2t{]:?xl02° y. Using 
the phase-space factor value G2!/=1.64 x 10~^^ y~^ (for g^i = 1.254) and the 
measured half-life, one obtains the nuclear matrix element value for the 21/(3(3 
transition M2v(0Ji'")=0.105 ± 0.010 (scaled by the electron rest mass). 

Because the technique used in [19] and in the present work does not allow 
a distinction between Ou^P and 2i//3/3 decay, strictly speaking, our result for 
Ti/2 is the sum of the 01/(3(3 and 2t//3/3 decay half-life times. However, the 
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experimental limit for the Oi>(3(3 decay to the O"*" ground state of ^"^Ru is 
about three orders of magnitude larger [20] than the value reported here for 
'^i/f'^^^''^^- Therefore, considering the reduced phase space available for the 
transition to the excited Of state, it is safe to assume that our result for T1/2 
refers solely to the 2z//3/3 decay. Additionally, theoretical estimates for Oiy(3(3 
decay to this level for a neutrino mass of 1 eV are many orders of magnitude 
higher: (7.6 - 14.6) x lO^^ y and 2.6 x lO^^ y [27]. 

In Ref. [15 the half-life ratio, Rt = Ti/2{0f) /Ti/2{0'^,s.) calculated for 
lo^Mo using the single-state dominance approach of [2B] and applying differ- 
ent scenarios for neutrino statistics. For purely fermionic neutrinos this ratio 
turned out to be Rt = 61 while for purely bosonic neutrinos the value Rt = 73 
was calculated. Using the very precise result of Ti/2=(7.11±0.54) xlO^^ y from 
the NEMO-3 experiment [20] for the O"*" ground-state transition and our half- 
life for the of excited state transition one obtains the ratio Rt = 77lig. Al- 
though this ratio seems to favor the hypothesis of bosonic neutrinos, a more 
accurate half-life time for the Oi" transition is needed to draw any definite 
conclusions. 



4 NEW DATA FOR DOUBLE-BETA DECAY OF ^^^Mo TO HIGHER 
EXCITED FINAL STATES OF i°0Ru 



Following [T9] we also determined new limits for the half lives of decays to 
higher excited states in ^°°Ru (see Fig. [5]). For the 2z//5/5 decay transition to 
the second excited 2+ state in ^°°Ru (22^) at 1362.2 keV we searched for the 
coincident detection of 822 keV and 540 keV 7 rays. No events of this type 
were observed. 

For the 2i//?/? decay transition to the second excited 0+ state of ^°°Ru (O2") at 
1741 keV we did not observe any events for the 1201 keV - 540 keV coincidence 
pair. For the 379 keV - 1362 keV coincidence pair one event was observed. 
This number is consistent with the background yield. 

In the case of the transition to the third excited 0+ state (Of) at 2051 keV our 
data show one 1512 keV - 540 keV event and no 689 keV - 1362 keV events. 
The latter event is consistent with the background yield. 

Finally, for the transition to the fourth excited O"*" state in ^°'^Ru (Of) at 2387 
keV we did not observe any events for the 1848 keV - 540 keV or the 1025 
keV-1362 keV coincidence pairs. 

Using the measured coincidence efficiency of Sec. 2 and correcting it for the en- 
ergy dependent attenuation in the ^°°Mo sample and energy dependent HPGe 
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detector efficiencies, new lower limits for the half-life time of the 2v(3(3 decay 
transitions discussed above were obtained. We combined these results with 
those of [19j using our newly determined coincidence efficiency and obtained 
the limits summarized in Tabled Basically, the limits were improved by about 
a factor of two compared to the earlier work [19] . 

The new limit for the 22 state transition is of special importance regarding 
the issue of neutrino statistics. Contrary to the 0^" state discussed in Sec. 3, 
the ratio Rt calculated using excited 2"*" states instead of excited 0"*" states is 
very different for bosonic and fermionic neutrinos. According to [H] the ratio 
Rt decreases for bosonic neutrinos by two orders of magnitude for the ffist 
excited 2^. Of course, this state is inaccessible with our coincidence technique, 
but the 22 state certainly is, and it is conceivable that our present limit of 
Ti/2 >4.4xlO^^ y for this state can be improved considerably in a dedicated 
experiment. 



5 SUMMARY AND CONCLUSION 



We performed an improved measurement of the coincidence detection effi- 
ciency of the TUNL-ITEP double-beta decay apparatus. The efficiency has 
been determined to an overall accuracy of 5.1%, providing the basis for deter- 
mining 2z//3/3 transitions to excited states in daughter nuclei to higher precision 
than previously possible. We also present new data obtained during a 450 day 
measurement of T1/2 for the 2v(3(3 decay of ^^''Mo to the ffist excited 0^^ state 
in ^°°Ru. Combining this data set with the previous data of [19] resulting in 
a total accumulation time of 905 days, we obtained the most accurate value 
for the half-life and NME of this decay. This half-life and our lower bound 
for the transition to the 22 state in ^°°Ru could be used to set bounds on a 
speculated bosonic component to neutrino statistics. 
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Fig. 1. Pulse height of HPGe detector 2 versus pulse height of HPGe detector 1 
with the events associated with the 468.6 keV - 475.1 keV decay sequence of ^'^^Ru 
shaded in. 

[26] F. Simkovic et al., Phys. Rev. C 64 (2001) 035501. 
[27] J. Suhonen, Nucl. Phys. A 700 (2002) 649. 
[28] P. Domin et al., Nucl. Phys. A 753 (2005) 337. 
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Fig. 2. An example of a projection of the pulse height distribution from Fig. [T]onto 
the HPGe detector 2 axis and fit to data. 



Table 1 

Summary of systematic error contributions. 
Uncertainty Contribution % 
Intensity of Calibration Gamma Source 3% 
Energy and Attenuation Correction Factor 1% 
z-dependence Correction Factor 1% 
Geometry of -'^'^'^Mo source 3% 

Non-Symmetrical Efficiency Curve 2.4% 

Dead Time 0.15% 

Uncertainty in ^^^Rh half-hfe 0.15% 

Total 5.1% 
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Fig. 3. Average of coincidence detection efficiency data obtained for the horizontal 
and vertical scans as a function of r for the E^i=590.8 keV and E^2=539.5 keV coin- 
cidence. The data were taken with the^^^Rh source (E^i=468.6 keV and E^2=475.1 
keV) and then corrected for detection efficiency and attenuation differences between 
the two 7-ray pairs involved. The curve through the data presents a least-squares 
fit. The upper and lower curves indicate the ±5.1% scale uncertainty associated 
with our data. 



Table 2 

Summary of half-life times and limits for the 2vj3l3 decay of ^'^'^Mo to various excited 
states of ^'^'^Ru. Limits are given at the 90% C.L. The results are based on the 
combination of the previous data [19J and the present data. See Table 11 of Hornish 
et al. 1191 



Transition Level (keV) Qpp (keV) T^J+^i- (x 10^° yr) T^J+^i- (x 10^° yr) from [E 

5.5lo:8 (stat) ± 0.3 (syst) 6.0l{:? (stat) ± 0.6 (syst) 
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Fig. 5. Level scheme of Ru. Only levels and transitions of interest to the present 
work are shown. 
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Fig. 6. These spectra show our recent cohected events over 450 days, (a) 7-ray 
spectrum in coincidence with the 590.8 it 3.0 keV transition in ^'''^Ru, and (b) 
7-ray spectrum in coincidence with the 539.5 it 3.0 keV transition in ^'^'^Ru. 
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Fig. 7. These spectra show the total collected events with the TUNL-ITEP appa- 
ratus over 905 days, (a) 7-ray spectrum in coincidence with the 590.8 =b 3.0 keV 
transition in ^'^'^Ru, and (b) 7-ray spectrum in coincidence with the 539.5 it 3.0 
keV transition in ^'''^Ru. 
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